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ABSTRACT 

We construct a two-zone model to describe hydrogen and helium burning on the surface of an accreting neu- 
tron star and use it to study the triggering of type 1 X-ray bursts. Although highly simplified, the model repro- 
duces all of the bursting regimes seen in the more complete global linear stability analysis of Naravan & Hevll 
J2003I) . including the regime of delayed mixed bursts. The results are also consistent with observations of 
type I X-ray bursts. At low accretion rates M/MeiIA ^ 0.1, thermonuclear helium burning via the well-known 
thin-shell thermal instability triggers bursts. As M increases, however, the trigger mechanism evolves from 
the fast thermal instability to a slowly growing overstability involving both hydrogen and helium burning. The 
competition between nuclear heating via the /9-limited CNO cycle and the triple-a process on the one hand, and 
radiative cooling via photon diffusion and emission on the other hand, drives oscillations with a period approx- 
imately equal to the hydrogen-burning timescale. If these oscillations grow, the gradually rising temperature 
at the base of the helium layer eventually provokes a thin-shell thermal instability and hence a delayed mixed 
burst. This overstability closely resembles the delayed mixed bursts of Narayan & Heyl. For M/Mem }t 0-25 
there is no instability or overstability, and there are no bursts. Nearly all other theoretical models, including de- 
tailed time-dependent multi-zone calculations, predict that bursts should occur for all M/MEdd ^ 1^ in conflict 
with both our results and observations. We suggest that this discrepancy arises from the assumed strength of 
the hot CNO cycle breakout reaction '^0(a,7)''^Ne in these other models. That observations agree much better 
with the results of Narayan & Heyl and our two-zone model, both of which neglect breakout reactions, may 
imply that the true '^0(a,7)''^Ne cross section is much smaller than assumed in previous investigations. 
Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances — stars; neutron — X-rays: 
binaries — X-rays: bursts 



1. INTRODUCTION 

Type I X-ray bursts are thermonuclear explosions that occur 
on th e s urfaces of accretin g neutro n stars (Babushkina et aQ 



T97l iGrindlav&Heisd fT97l Icirindlav et al. 
BelianetalJ Il976t IWooslev & TaanJ 119761: LTossI [197^ 

Maraschi & Cavaliere '1977: 'Lamb & Lamb' '1977. 



They are triggered by unstable helium or hydrogen burning 
(for reviews, se e Lewin et al. 1993, 1995; Cumming 2004; 
IStrohmaver & Bildstea.2006) . For a wide range of accretion 
rates, the basic physics of the burst onset is that of the well- 
known thin-shell thermal instability ('Sc hwarzschild & HarmI 
[1965; Hansen & van Horn 1975), which operates as follows. 
In equilibrium, the energy that is generated by nuclear 
burning is exactly balanced by the outward diffusion of 
photons coupled with radiation from the stellar surface. If the 
nuclear reaction rates are sufficiently temperature-sensitive, a 
positive temperature perturbation causes the nuclear energy 
generation rate to increase faster than the rate at which the 
stellar surface can cool, causing a thermonuclear runaway 
and hence a type I X-ray burst. The physics of this instability 
can be described quite well with a simple one-zone model 
jpuiimoto et al. 1981; Paczvnski 1983a; Bildsten 1998; 
ICumming & Bildsten 2000). If the amount of accreted 
matter needed to ignite a burst were roughly independent of 
the accretion rate M, as suggested by simple models, one 
would expect the recurrence time fiec to vary as M"'. This 
is in reasonable agreement with observations for a range of 
M < O.lMfidd, where Mndd denotes the mass accretion rate 
at which the accretion luminosity is equal to the Eddington 
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limit. However, for 0.1 < M/MEdd 0.3, observations show 
that free increases with increasing M, and bursts are found 
to ce ase altoge ther for M /MpM ^ 0.3 ( van Paradiis et aO 
119791 fl988; Cornehsse et al. .2003; .Remillard et al. 2006I) . 
Moreover, for 0.1 < M/MEdd ^0.3, most of the accreted 
matter appears to burn stably between consecutive bursts 
(lyanParadijs et al. 1988; in't Zand et al. 2003). Simple 
thin-shell thermal instability models, on the other hand, 
predict that bursts should occur for all M/Medd ^ 1 and that 
tr er. should dccreasc mono tonically with increasing M. 

iNaravan & Hevl (2003, hereafter NH03) developed a 
global linear stability analysis of the accreted nuclear fuel 
on the surface of a neutron star. Unlike thin-shell thermal 
instability models, for which one assumes that purely ther- 
mal perturbations trigger bursts, NH03's analysis was quite 
general and made no assumptions as to the trigger mech- 
anism. Apart from confirming the standard thermal insta- 
bility for M/Medd ^0.1, NH03 discovered a new regime 
of unstable nuclear burning for 0.1 < M/MEdd ^0.3 that 
they referred to as "delayed mixed bursts." This regime re- 
produced the increase in free with increasing M seen in ob- 
servations, as well as the occurrence of considerable sta- 
ble burning between successive bursts. It also predicted the 
absence of bursts for M/Mp ad ^, 0.3, again in agreement 
with o bserva t ions. Recently, Cooper & Naravaa ( i200 5) and 
Coop er et alJ J2006h have suggested that the stable burning 
that precedes delayed mixed bursts is crucial for the produc- 
tion of the carbon fuel ne eded to trigger superbursts (for re - 
views of superbursts, see Kuulkers 2004; Cumming 2005^, 
since only stable nuclear burning can produce carbon in suf- 
ficient quantities JSchatz et alJl999ll200!ll2003tlKoike et alJ 
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l2()04llWooslev et alJ2Q04tlFisker et alJ2005h . 

It is perhaps not surprising that a global linear stability anal- 
ysis is able to reproduce type I X-ray burst observations better 
than simple thin-shell thermal instability models. However, 
the fact that the results of all time-dependent multi-zone mod- 
els with large reaction networks c onflict with observations 
as wen (AvasH&Joss 1982; Taa m et alJll99d: iFisker et al. 
l2()03trHegeret al..2005.) is much more problematic. Like the 
global linear stability analysis, these latter models include no 
preconceived instability criteria. Therefore, they should in 
principle be able to reproduce the regime of delayed mixed 
bursts and the absence of bursts for M/MEdd ^ 0.3, as ob- 
served in nature. The disagreement between the predictions of 
time-dependent multi-zone burst models and the global linear 
stability analysis of NH03 suggests that some crucial aspect 
of the physics is different between the two. Understanding 
this discrepancy is of fundamental importance in the study of 
type I X-ray bursts. 

While the global linear stability analysis of NH03, with its 
identification of the delayed mixed burst regime, succeeded in 
reproducing several puzzling observations of bursts, unfortu- 
nately, the complexity of the model and the absence of a spe- 
cific stability criterion make it difficult to understand the basic 
physics of delayed mixed bursts. To rectify this situation, we 
construct here a simple two-zone model that reproduces all of 
the basic features of the global analysis. We begin in |2]with 
a derivation of our model. In ^ we discuss the equilibria 
of the model and analyze their stability as a function of ac- 
cretion rate. We then integrate the equations that govern the 
model and use the results in l^to explain the physics of the 
onset of delayed mixed bursts. We discuss our results in ^Sj 
and we conclude in ^ 

2. THE MODEL 

We assume that matter accretes spherically onto a neu- 
tron star of gravitational mass M = 1 AMq and areal radius 
R = 10 km at an accretion rate per unit area E as measured in 
the local frame of the accreted plasma. We consider all phys- 
ical quantities to be functions of the column depth S, which 
we define as the rest mass of the accreted matter as measured 
from the stellar surface divided by AttR^. We denote the Eu- 
lerian time and spatial derivatives as d/dt and d/dY., respec- 
tively, and we denote the Lagrangian derivative following a 
parcel of matter as D/Dt, where D/Dt = d/dt + td/dY.. We 
describe the composition of the matter by the hydrogen mass 
fraction X, helium mass fraction Y, and heavy element frac- 
tion Z = l-X-Y. The mass fractions at S = 0, Xq, Yq, and 
Zo, are those of the accreted plasma. In this work, we assume 
that all heavy elements are CNO and that the composition of 
the accreted matter is that of the Sun: Xq = 0.7, Yq = 0.28, 
Zo = 0.02. 

The following set of four partial differential equations gov- 
erns the thermal structure and time evolution of the material 
(e.g., NH03; see Table 1 of that paper for the definitions of 
the various symbols): 
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Since we are interested in developi ng a sim ple model that 
avoids unnecessary details, we follow lBildsten, U.99^ and set 
the opacity to a constant value, k = 0.4Kes = 0.136 cm^g"', 
where Kes is the Thomson scattering opacity. We assume that 
a nondegenerate, ideal gas supplies the pressure P, which is 
a good approximation for the shallow column depths we con- 
sider in this work. The specific entropy s is in general a func- 
tion of both T and P = gY, and so 



Ds 
T — : 
Dt 



dt 



dY E 



(5) 



where we set the specific heat at constant pressure Cp = 
5ks/2mp for an assumed mean molecular weight /i = 1, and 
Vad = {dlnT / dlnY)s = 2/5 for an ideal gas. We ignore the 
terms oc S for reasons we will give in ^(this is a standard ap- 
pro ximation in most one-zone models of type I X-ray bursts, 
e.g. jBildstenlll998tlCumming & Bildstenl2000h although we 
have done calculations with these terms included, and we have 
verified that they are not important. Thus, equation (O be- 
comes 

dF dT 

^=Cp — -(eH+EHe). (6) 



dt 



__The usual one-zone approach to type I X-ray bursts 
(Fuiimotoetal. 1981; Paczvhskil Il983at iBildsterJ 119981; 
■Gumming & Bildsten 2000) focuses primarily on helium 
burning, for which a single zone is adequate. Several authors 
have introduced more elaborate two-zone models of type I 
X-ray bursts (Buchler & Perdang 1979; Barranco et al. 198(J 
IRegev&Livio 1984; Livio & Regev 1985; Yasutomi 1987), 
but these models usually include only helium burning, and all 
of the nuclear burning is confined to a single zone. Jos e et alJ 
i 19931) present a two-zone model of thermonuclear burning on 
accreting white dwarfs that includes both hydrogen and he- 
lium burning, but the two fuels are forced to burn in separate 
zones by construction. This approximation is too severe for 
realistic modeling of type I X-ray bursts. The global linear 
stability analysis of NH03 clearly shows that the interaction 
between hydrogen and helium burning plays a key role in de- 
layed mixed bursts. Motivated by this, we describe in this 
paper a new model with the following two zones: (i) a zone 
that begins at the surface of the star, where S = 0, and extends 
to the depth Eh at which hydrogen is depleted via nuclear 
burning, and (ii) a zone that begins at Eh and extends to the 
depth Ehc at which helium is depleted via nuclear burning. 
Hydrogen, helium, and CNO are all present in zone (i), while 
only helium and CNO are present in zone (ii). Furthermore, 
we include both hydrogen and helium burning in zone (i), but 
we include only helium burning in zone (ii) since hydrogen 
is not present. Note that, in our model, hydrogen fully de- 
pletes before all of the helium burns by construction, so that 
Ehc > Eh. This condition is satisfied over the entire range 
of E of interest in this work. We discuss the case for which 
Ehc < Eh in ^ We make the simplification that all ther- 
monuclear processes within a zone take place at the bottom of 
that zone. Equation Q thus implies that the flux F in a par- 
ticular zone is a constant with respect to E, and so T"* varies 
linearly with E within each zone (see equation^. 

In the simplest one-zone models, all parameters of the mat- 
ter are taken to be constant within the zone. Clearly this pre- 
scription is a gross simplification, for in reality the physical 
parameters X, Y, Z, T, and F vary smoothly as functions of 
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S. In order to improve the accuracy of our model, we assume 
linear profiles with respect to depth of the various quantities 
of interest (except F). Thus we take 



X(S],f) = 



^o(l- 
0, 



S>I]h, 



(7) 



Z(E,f): 



Zo + (Zh-Zo)(S/Sh), S<Sh 

Zh + (1-Zh)(E-Sh)/AS, I]H<S<SHe 
1, S>SHe, 

(8) 

Eh < S < Ehb 



7h'(S/Eh), 

r4+(r4^-r4)(s-SH)/AE, 



7^4 



E > E 



He, 




E<Eh 

Eh < E < Ene 

S > Ehe, 



(9) 
(10) 



where 



Zh=z(Eh), rH = r(EH), rHe = r(EHe), (11) 

AE = EHe-EH, (12) 

and the five quantities Eh, Zh, 7h, Ehb, and Thc are functions 
of time. Note that equation (|9j is exact for the aforementioned 
simplifications of constant opacity and constant flux within a 
zone. Figure^shows a pictorial representation of the model. 
For simplicity, we assume that the flux F entering the bottom 
of the helium-burning zone from the stellar interior is zero. 
Integrating equation Q gives 



Fh = 



acT^ 

3kEh ' 



Fup, — 



ac(Tl-T^) 



(13) 



(14) 



3kAE 

For the relatively high accretion rates we consider in this 
work, hydrogen burns to helium via the hot CNO cycle, so the 
corresponding nuclear energy generation rate = E^jnZ is a 
function only of Z, where 7h = 9. 1 x 10"'* s"' ( Hovle & Fowled 
[T965). Helium burns to carbon via the triple-a reaction. The 
corresponding energy generation rate ene is a strong function 
of temperature, but it depends also on Y and the density p. In 
this work, we simplify matters by approximating ene to be a 
function only of temperature (Hansen et al. 2004), so we set 

eHe(r) = £;^,7He(7i)/r)3exp(-ro/r), (15) 

where 7b = 4.4027 x lO^K and we choose 7He = 0.1 s~' such 
that EHe in our model is a reasonable representation of the true 
eiie(T, Y, p) at densities typical of those at which helium burns 
on accreting neutron stars. 

We wish to derive a set of five differential equations that 
determine the behavior of the five physical quantities Eh, Zh, 
Til, Ene, and The. To do this, we integrate equations Q, Q, 
and from to Eh and equations (0} and from Eh to 
Ehc- This gives 



dX ■ dX 

dZ -dZ 
dt aE 



JE = - 



di: = 



dT 



-^i^E, 



t/E, 



^He 



(16) 



(17) 



(18) 
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He 







-He 



Fig. 1. — A schematic representation of how the composition, temperature, 
and flux vary as a function of column depth in the proposed two-zone model. 
Xq, Yq, and Zq denote the hydrogen, helium, and CNO mass fractions of the 
accreted plasma, respectively. Zone (i) extends from E = to S = Sh, and 
zone (ii) extends from E = Sh to S = Ene- 



dt dY. 



dF 

aE^^ = 



dt ' 



EHe dT.. 



(19) 



(20) 



Since ene is a very strong function of temperature, most of 
the energy generated via helium burning within a zone will be 
released near the bottom of that zone, where the temperature 
is greatest. Therefore, when we integrate over ene, we make 
the approximation that the temperature throughout the zone is 
the temperature at the base, so we set 

enet/E = eHe(7H)SH, (21) 

/o 

eHe^^E EE eHe(7He)AE. (22) 

In contrast, en is a weak function of depth, and so the energy 
generated via hydrogen burning within zone (i) will be re- 
leased more evenly throughout the zone. Therefore, when we 



integrate over en, we perform the integral exactly according 
to equation (|8j: 

f^H /'Sh /y \ 

(23) 

To derive the system of differential equations that go vern our 
model, we evaluate the integral s in eq uations ( I16I2Q> accord- 
ing to the ansatze of equations J7I10> and the approximations 
described above. Using the ansatz forX(I],f) in equation Q 
and the result of equation ( I23t . equation il6\ becomes 



di 



^0 1 





■7h 



as 

Zh + Z() 



(24) 

From the fundamental theorem of calculus, and recalUng that 
Sh = SH(f)i this gives 

i-.— . ^^^^ 

By performing analogous steps on equations ( I17I20K equa- 
tions ( I16I20> may now be written as 



1 ^Sh ■ /Zh + Z{) 



(26) 



— Sh-(Zh-Zo) — 



eHe(7H) 



Sh-(Zh-Zo)S, 



^He 



(27) 



€He(7H)SH - (fk - ^He) , 

(28) 



/asHe aEe 

— AS-(1-Zh)( ^ + ^ 



eHe(7He) 



AI]-(1-Zh)S, 



(29) 



"He 



C„ I \ ' V l] = eHe(7ke)AS-fiie. (30) 



9f 



As the quantities Sh, Zh, 7h, She, and Tee are func tions of 
only one independent variable t, equations ( I26I30> can be 
rewritten as the following set of coupled ordinary differential 
equations: 



^^«.=2S-^(Zh + Zo), 



dt 



dZH ^ £He(7H) _ (Zh-Zq) / ^ t/Ee 



dt 



dt 4Cr, 



"He 



'E'h7H ( ^' I +eHe(7H) 



dt 

7h ^Sh 



(31) 
(32) 



(33) 



He 



21]- 



(1-Zh) 



, eHe(7He) dZyi 

dt~ 



"He 



(34) 



dT^ 



He 



5(7H'e- 



4\2 



7h^) 



dt 



4C,(Tl-5T^X + 4TlTi) 



™ . -pHe Cp(Tiis-Tii) d'En 



^He 

"4AE 



Tke 



c/AS 
dt 

dTii 
dt 
(35) 



These are the fundamental equations of our two-zone model. 
They can be expressed in the compact form 



dx 
dt 



= f 



S,Xo,Zo 



(X), 



(36) 



where x represents the vector {Sh,Zh, 7h, ^He, Tne}, and S, 
Xq, and Zo are free parameters. 

The reader should note that several of the dimensionless 
coefficients in the above equations are ~ 1 and may be set 
to unity. Furthermore, Eh and Ene change over much longer 
timescales than Tu and Tne, and so the d'Eu/dt and d'Sjie/dt 
terms in equations ( I33t and i35\ can be set to zero. Therefore, 
equations i33i and i35i can be replaced by 



dTn _ 
dt ' 4C 



j^* , Zh+Zo\ ^H--pHe 

^h7h 1 — — ) +eHe(7H) 



p 

dTn, _ 5 

dt 4Cn 



eHe(7He)" 



AE 



Sh 

dTn 

' dt ' 



, (37) 



(38) 



with little loss of accu racy, al though the results we present are 
for the full equations ( I31I35> . 

3. EQUILIBRIA AND THEIR STABILITY 

Our goal in this section is to find the equilibrium solutions 
of the above two-zone model and to determine their stabil- 
ity. To fin d the eq uilibrium solutions, we set d/dt = 0, so that 
equations ( I31I35> become the following set of 5 coupled al- 
gebraic equations: 



7h 



Zh + Z(i 



EXo 



eHe(7H)EH = E(ZH-Zo)£'He; 
FH-FHe = £;57H ( ) EH + eHe(7k)SH, 



eHe(7ke)AE=E(l-ZH)£He, 
Fue = eHe(7He)AE. 



(39) 

(40) 

(41) 

(42) 
(43) 



We solve this set of equations for the 5 equilibrium values 
Ey , Z'^, T^'^, E^j., and T^^. When we do this, we obtain 



F^'f = E[Xo£^ + (l-Zo)£ii,], 

F^^ = E(l-Z^'')£He- 



(44) 
(45) 



Thus the flux emitted from the stellar surface, F^'*, equals the 
flux released via steady-state nuclear burning of all the fuel, 
and the flux entering zone (i), F^, equals the flux released 
via steady-state nuclear burning of the helium within zone 
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log(Lec) 



Fig. 2. — Equilibrium values of the five variables Th, Tho Sh, Sho and Zh 
as a function of /acc = S/SEjd- Note that Tj^^' > T^^ and > Ejf' for the 
indicated range of /acc ■ For larger accretion rates, the inequalities are reversed 
and the model described in this paper is no longer valid. 



(ii). Note that the equilibria are functions of E, Xq, and Zq. 
See Figure 121 for plots of the equilibrium values of the five 
fundamental variables as a function of the Eddington-scaled 
accretion rate 

t 



/aec=-^, (46) 
^Edd 

where Endd = 1.0 x 10^ gcm"^s~'. 

The behavior of the equilibrium quantities shown in Fig- 
ure |2]is easily understood. Nuclear burning generates energy 
and heats the accreted layer. In equilibrium, hydrogen and he- 
lium burn at a rate oc 1.^^^, and so the equilibrium temperatures 
T^'^ and TJ^^ are monotonically increasing functions of /acc- 
The rate at which helium burns is very temperature-sensitive, 
and so the depth at which helium depletes, E^^., is a decreas- 
ing function of U^c- One-zone models often presume that no 
stable helium burning occurs prior to a burst, in which case 
= Zfi. This wou ld imply that E^ oc Ig^^c from equation 
J39> . lTaam & Picklum (1978) were the first to remark that sta- 
ble helium burning prior to ignition increases the abundance 
of seed nuclei for the hot CNO cycle and thus expedites hy- 
drogen burning. Figure |2l shows that this effect is significant 
for the accretion rates we consider, since is often much 
greater than Zq. Consequently, far from increasing with /acc, 
E^'' is in fact a decreasing function of /acc- 



> 
a 

W 



> 



100 
80 
60 
40 
20 

-20 
-40 

1 

0.8 
0.6 
0.4 
0.2 






^('^max^'-acc 
'^H tacc 
















Th 











-1.4 -1.2 



-1 -0.8 

log(lacc) 



-0.6 -0.4 



Fig. 3. — Top panel: spectrum of the real parts of the eigenvalues as a 
function of Zacc- Also shown are the imaginary part of the eigenvalue with 
the greatest real part and ujn = In/tn, where fn = S^'/S is the hydrogen- 
burning timescale. Each frequency is multiplied by the accretion timescale 
?acc = S^1./E. Bottom panel: normalized components of the eigenvector 
corresponding to the eigenvalue with the greatest real part. 



After we have solved for the equilibrium solution, we 
conduct a linear stability analysis to determine whether 
the system is stable or unstable to perturbations (e.g., 
Rjuckenheimer & Holmes 1983). From equation i36l . 



d(Sx) 
dt 



(47) 



where dt/dx denotes the Jacobian matrix, which is evalu- 
ated at the equilibrium x'"*. We calculate the eigenvalues 
of the Jacobian numerically us ing the standard procedure 
for a real, nonsymmetric matrix (Wilkinson & Reinscy il97ll 
Press et al. 1992). If at least one of the eigenvalues has a posi- 
tive real part, then the system is unstable to perturbations, and 
we say that type I X-ray bursts occur. Otherwise, the system 
is stable to bursts. Thus the eigenvalue with the greatest real 
part determines the bursting behavior of the system. 

The top panel of Figure |3l shows the spectrum of the real 
parts of the eigenvalues A as a function of /acc- We have mul- 
tiplied the eigenvalues by the accretion timescale. 



(48) 



to make them dimensionless. The most negative eigenvalue 
has 5R(A)facc S -1000 and is not shown. The bottom panel of 
Figure |3l shows the normalized squared moduli of the com- 
ponents of the eigenvector v corresponding to the eigenvalue 
with the greatest real part. To construct v, we nondimension- 
alize the perturbations 6xi such that the components of the 
eigenvector are v,- = (51n(x,). 

Figure nicely illustrates the nature of the instability and 
the behavior of the resulting type I X-ray burst as a func- 
tion of the accretion rate. For log(/acc) S -1-3, A is real and 
V is dominated by perturbations in Thc, showing that a pure 



6 



thermal perturbation in the heUum burning region triggers the 
instability. Also, the perturbation grows very quickly since 
^Wt-dcc ^ 1- Hydrogen burning is clearly inconsequential to 
the instability since perturbations grow at a rate K(A) > ojh, 
where 

c^H = — (49) 

and 

v^eq 

tH = ^ (50) 

is the hydrogen-burning timescale. These bursts (named he- 
lium bursts by NH03) are triggered by the well-known and 
standard thin-shell thermal instability. 

For -1.3 < log(/acc) < -0.9, 5R(A)4cc > 1 so the bursts 
are still "prompt," but now A is complex. This means that 
there is an overstability rather than a pure instability. Also, 
5R(A) < LUii, which means that hydrogen burning now begins 
to play a role in the stability of the system. As a confirmation, 
we note that v involves perturbations of all five parameters, 
in particular the composition of the hydrogen-burning zone, 
showing that all the variables play a role in the instability. 

For -0.9 < log(4cc) < -0.6, K(A)facc < 1 and 5(A)/K(A) > 
1. This is the regime of delayed mixed bursts identified by 
NH03. The growing mode is overstable and undergoes sev- 
eral oscillations as the amplitude slowly increases. The an- 
gular frequency of the oscillations 5(A) « luh, showing that 
hydrogen burning plays a key role. We have verified that the 
imaginary parts of the eigenvalues in this regime derived us- 
ing the global linear stability analysis of NH03 are approxi- 
mately equal to cjh as well. 

Finally, all eigenvalues have negative real parts for 
log(/acc) ^ -0.6, and so no bursts occur at these high ac- 
cretion rates. The See component begins to dominate v in 
this regime. For these systems, the effective radiative cool- 
ing rat e of the accreted la yer, which is proportional to 
(e.g., Fuiimoto et al. 1981; Bildsten 1998), is large enough to 
dampen any oscillations that might otherwise initiate an insta- 
bility. 

It is now clear why we are able to drop the terms oc E in 
equation (|2}- In the helium and prompt mixed burst regimes, 
these terms are unimportant because thermal perturbations 
drive the instability, and they occur on a timescale much 
shorter than the accretion timescale t^^c (e.g., Bildsten 1998). 
In the delayed mixed burst regime, although the mode grows 
on a timescale comparable to the accretion timescale, compo- 
sitional perturbations primarily drive the overstability. 

We consider our five-parameter, two-zone model to be 
"minimal," by which we mean that no model consisting of 
only a proper subset of the five time-dependent variables Eh, 
Zh, 7h, Ehb, and Tee is able to reproduce the phenomenon of 
delayed mixed bursts. To demonstrate this, and to illustrate 
the specific role each parameter plays, we conduct linear sta- 
bility analyses on truncated n-parameter models, where n < 5, 
and we compare the results of these models to those of the full 
model. To perform these calculations, we determine the equi- 
libria of the five parameters as usual, but we perturb only the n 
parameters Sxj and set the perturbations of the other 5 - « pa- 
rameters 6xi = 0. We then solve for the eigenvalues of the cor- 
responding nxn Jacobian matrix. See Figure0]for plots of the 
real part of the largest eigenvalue of each model as a function 
of Zacc- The middle dotted line is the eigenvalue for an isobaric 
purely thermal perturbation, where (5Eh = ^Ehb = SZu = 0. All 
equilibria of this model are stable for log(/acc) ^ -1.2. This 




-1.4 -1.2 -1 -0.8 -0.6 

log(lacc) 

Fig. 4. — Real part of the largest eigenvalue as a function of tec for three 
models in which only a proper subset of the five parameters Eh, Zh, 7h, Shc, 
and 7hc is perturbed. The dotted lines show the normalized eigenvalues for 
models in which from top to bottom (i) (5Sh = <5Shc = 0, (ii) (5Sh = <5Shc = 
(5Zh = 0, and (iii) (5Zh = 0. The solid line shows the normalized eigenvalue 
for the full model. 



result nicely complements the eigenvector analysis shown in 
Figure|3l which illustrates that thermal perturbations drive the 
instability for log(4cc) ^ -1 -2, but they become less important 
for higher accretion rates. Using their global linear stability 
analysis, NH03 performed similar calculations in which they 
considered only thermal perturbations, and they obtained very 
similar results (see their Figures 16 and 17). The bottom- 
most dotted line is the eigenvalue for a model in which only 
(5Zh = 0. This is a thermal perturbation as well, i.e., it includes 
no perturbation in the nuclear composition parameter Zh, but 
it is one in which the isobaric constraint is relaxed. Figure |4] 
shows again that all equilibria are stable for log(Zacc) S -1-2. 
However, the eigenvalues near the transition between stabil- 
ity and instability are now complex. This is evident from the 
kink in the line at log(Zacc) ~ -1.35, indicating that two real 
eigenvalues have merged into a complex-conjugate pair. The 
topmost dotted line is the eigenvalue for a model in which 
(5Eh = (5EHe = 0. This perturbation is again isobaric, but it in- 
volves both thermal and compositional perturbations. Figure 
13 shows that the eigenvalue is real, and that all equilibria are 
unstable. 

To better understand the role each parameter plays in the 
stability of the dynamical system, we first divide the perturba- 
tions into two classes. Equations ( I33> and i35i imply that STu 
and STus change on timescales on the order of the thermal dif- 
fusion timescale, whereas equations ( 13 U . i34l . and i32\ imply 
that (5Eh, (JEhb, and dZu change on timescales on the order of 
the accretion timescales tu and t^^c, which are usually much 
longer than the thermal diffusion timescale. Figures |3] and |4] 
demonstrate that, for log(/acc) ^ -1-2, purely thermal pertur- 
bations drive the instability, which grows on a timescale much 
shorter than the accretion timescale. A positive temperature 
perturbations STu^ > grows quickly and accelerates helium 
burning, which decreases (JEhb- Although a decrease in (JEhc 
increases the effective cooling rate, it does so over a much 
longer timescale, and so the effective cooling rate is usually 
unable to compensate for the positive temperature perturba- 
tion. Therefore, the system is still thermally unstable. Thus, 
for thermal perturbations, it makes little difference whether 
or not one allows (5Eh and (JEhb to vary with regard to the 



Two-Zone Burst Model 



7 



stability of the system (l Paczvnskil983aHHeger et al.l2005l) . 

The situation is quite different for the delayed mixed burst 
regime, in which compositional perturbations 6Zu primarily 
drive the instability. A positive perturbation of any one of the 
three parameters 6Tu, STue, or SZu will cause the other two 
to grow as well, for a positive temperature perturbation accel- 
erates helium burning, which produces more CNO and hence 
increases 6Zu, and a positive 6Zu perturbation accelerates hy- 
drogen burning, which raises the temperature and hence in- 
creases 5Tu and STus. Since this system is stable to purely 
thermal perturbations, any unstable modes grow on the slow 
accretion timescale. But JSh and JEhc, which regulate the 
effective cooling rate, vary over a similar timescale. There- 
fore, JSh and 6J^ue have time to respond to changes in the 
nuclear heating rate in this regime. A positive perturbation in 
any one of 6Tu, STue, or SZu grows slowly and accelerates 
both hydrogen and helium burning, which decrease SJ^h and 
(JEhc- Decreases in SJ^h and (JShe raise the effective cooling 
rate over a similar timescale, which lowers the temperature 
and dampens nuclear burning. Consequently, SJ^u and SJ^ue 
then slowly increase as freshly accreted matter advects into 
the layer until nuclear burning begins again. This cycle re- 
peats, driving oscillations with a period on the order of the 
hydrogen burning timescale fe- 

4. TIME EVOLUTION OF THE ONSET OF DELAYED MIXED 
BURSTS 

Although linear stability analyses of the set of governing 
differential equations clearly determine whether or not the 
system will exhibit type I X-ray bursts, they can only hint 
at the nonlinear development of the instability and the un- 
derlying physics of the burst o nset. T o remedy this, we 
numerically integrate equations ( I31I35> to study the non- 
linear evolution of the onset of delayed mixed bursts. We 
will not discuss the onset of prompt bursts in any de- 
tail, since the physics of prompt bursts has been previously 
studied by many authors and is generally well understood 
Cpss"T97l 'Taam&PicklumllT97P, 'Joss & Li' 198^, 
1980; AvasH & Joss 1982; Hanawa & Sugimoto 1982; Taa3 
1982; Paczvnski 1983bt iTaam et al.ll993tiZingale et al.»200U 
Woosley et al. 2004). 

We perform the integration for a system with log(/acc) = 
-0.7, which is well inside the delayed mixed burst regime. 
Figure |3] shows that the system is marginally unstable, with 
5?(A)facc ~ 0.75, and that the system undergoes oscillations 
with a period approximately equal to the hydrogen-burning 
timescale tu- We initiate the overstability by perturbing the 
equilibrium such that the initial conditions of the integration 
xo = l.OOlx'^1. Figure |5] shows the results of this calcula- 
tion. The top panel illustrates the time evolution of the phys- 
ical parameters. All modes other than the principal mode are 
transient and quickly decay after the initial perturbation. 7h 
and Thb are essentially in phase since the thermal diffusion 
timescale between Eh and Ehb is much smaller than fn- The 
other physical parameters are pairwise out of phase, but each 
slowl y grows and oscillates with the principal mode ( Buchler 
Il993h . One can ascertain the relative phases by contrasting 
the phases of the various components of the complex eigen- 
vector V. For a comparison, we plot the time evolution of the 
physical parameters for a sequence of prompt bursts from a 
system with log(Zacc) = -1-4 in Figure |6l The results of this 
prompt burst calculation are very similar to those of previous 

one- and two-zone models CBarranco et al 1980; iPaczvilskil 

fT983aHRegev & Livioll984tlHeger et al.l2005h . 




Fig. 5. — Example of the time evolution of the overstability preceding 
a delayed mixed burst. Xq = 0.7, Zq = 0.02, and log(/acc) = —0.7 for this 
calculation. Top panel: evolution of the normalized physical quantities. 
Middle panel: normalized lightcurve. Bottom panel: evolution of the ratio 
(d^Hc / dT) / (decoo\/ dT) evaluated at Ehc- A thin-shell thermonuclear insta- 
bility ensues when (deiic/dT)/(decoo[/dT) > 1. 



To better understand the physics of the overstability preced- 
ing a delayed mixed burst, we consider a positive pertur- 
bation. Since the thermal diffusion time from Eh to Ehc is 
very short, Tus follows the time dependence of Th through- 
out the integration. An increase in Tu causes the helium in 
zone (i) to burn at a higher rate, which increases Ze, and the 
simultaneous increase in The causes the helium in zone (ii) to 
burn at a higher rate also, which decreases Ehc- The larger 
Zh causes hydrogen to burn at a higher rate, which decreases 
Eh. This raises the effective cooling rate of the accreted layer, 
which is proportional to 1 /E^. When nuclear burning has 
nearly depleted the hydrogen in zone (i), both of the zones 
cool as freshly accreted matter advects down to larger col- 
umn depths. Tu and Tue consequently decrease until enough 
nuclear fuel has accumulated to restart nuclear burning. Once 
this happens, Tu and Tm will increase and the cycle continues. 

If the amplitudes of the oscillations increase in time, the 
temperature Tue at Ehc may raise the helium nuclear energy 
generation rate high enough to trigger a thin-shell thermal in- 
stability and hence a type I X-ray burst. To understand when 
the burst will occur, we follow the basic one-zone, single- 
parameter hnear stability analysis of iFuiimoto etail lfT98l 
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Fig. 6. — Example of the time evolution of prompt bursts. Xq = 0.7, 
Zq = 0.02, and log(/acc) = -1.4 for this calculation. The top panel shows the 
lightcurve, and the bottom three panels show the time evolution of the live 
physical quantities. Note that Tne > 7h and Shc > Sh. 



and consider the effective radiative cooling rate at Ehb 
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(d^He/ dT) /{decooi/ dT) in the bottom panel of Figure [S] 
A comparison between the lightcurve shown in the middle 
panel and the ratio (deHe/dT)/{d€cooi/dT) in the bottom 
panel shows that a short type I X-ray burst, characterized by 
a sudden increase in the outward flux, occurs when the thin- 
shell thermal instability criterion {deHe/dT)/(decooi/ dT) > 1 
is satisfied. Note that the low peak flux of the burst is due to 
the restriction in our model that hydrogen burns only via the 
temperature-independent hot CNO cycle. In our two-zone 
model, this delayed mixed burst calculation evolves to a 
limit cycle of short bursts with a recurrence time ~ t^. In 
rea lity, hydrogen burning wil l proceed via the rp-process 
of IWallace & WooslevI il981h during the burst, and this 



will burn nearly all of the fuel and produce an Eddi ngton- 
limited flux (e.g., p Tanawa & Fuiimoto 1984; Woosle v et alJ 
|2004|: iFisker et al.ll2006i) . The system will then restart near 



f/facc « in Figure|5]and undergo several oscillations before 
having its next burst. Our two-zone model, which includes no 
rp-process reactions, is too simple to reproduce this behavior 
Note also that, while equation (15 U is the exact cooling 
rate for a suitably constructed one-zone model, it is only 
approximate for our two-zone model. However, it is accurate 
enough to illustrate that, much like prompt bursts, helium 
burning ultimately induces a thin-shell thermal instability 
that triggers a delayed mixed burst. 

While the above calculation illustrates the physics of the 
onset of delayed mixed bursts, it does not unambiguously 
demonstrate that the bursts are in fact "delayed." Specifically, 
we have not yet demonstrated that a significant period of sta- 
ble burning will precede a type I X-ray burst. Although we 
began the integration by perturbing the equilibrium solution, a 
system that is dynamically unstable need not come anywhere 
close to its equilibrium. Unfortunately, we cannot set the ini- 
tial conditions to that of a bare neutron star, for which xo = 
{ Eh , Zh , 7h , Ehc , Tke } = { 0, Zq, 0, 0, } , because the time steps 
required to begin the integration from these initial conditions 
are too small to maintain numerical accuracy. Instead, we 



set Xo = 



7^=4 T-cq Y'^^ 



eq 



i.e. we start the system with 



the hydrogen-burning layer in place but without any helium- 
burning layer We then integrate until a type I X-ray burst 
occurs. This calculation gives a recurrence time free = tu+hnt, 
where t^t is the integration time. From the recurrence time, 
we derive an estimate for the dimensionless quantity a, which 
is defined as the accretion energy released between successive 
bursts divided by the nuclear energy released during a burst. 
Thus, 



[Xo£* + ( 1 - Zo)£iJe] + ( 1 - ZH)£j^eSHe ' 



(53) 



where z = (1 -2GM / Rc^)~^ 1^ - 1 is the gravitational redshift, 
and the values of Eh, Zh, and Ene are determined from the 
numerical integration just prior to the burst. Both the accre- 
tion energy and nuclear energy of the burst in equation ( 15 3> 
are those measured locally at the stellar surface. In our ex- 
pression for a, we implicitly assume that all of the hydrogen 
and helium burns to CNO during a burst. Figure shows 
that a ~ 100 in the regime of prompt bursts, which is consis- 
tent w ith more sophisticated calculations (e.g.. lWo oslev et al] 
l2004h . and that a rises dramatically to values ^ 1000 for 
log('acc) -0.9, which illustrates that the systems in the de- 
layed mixed burst regime undergo long periods of stable nu- 
clear burning prior to type I X-ray bursts. This is consistent 
with both the res ults of NH03 (see their Figures 15 and 18) 
and observations dvan Paradiis et alJl988h . 

5. DISCUSSION 

In this section, we compare and contrast the results of the 
two-zone model presented in this work with the results of both 
the global linear stability analysis of NH03 and other theoret- 
ical models. 

5.1. Comparison to the Global Stability Analysis of Narayan 

&Heyl 

The primary goal of this work is to construct a simple model 
of nuclear burning on accreting neutron stars that elucidates 
the essential physics behind the more complicated and ab- 
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Fig. 7. — Shows a-values of type I X-ray bursts as a function of /acc- The 
steep rise in q beginning at log(/acc) ~ -0.9 as /acc increases is characteristic 
of delayed mixed bursts. 



struse global linear stability analysis of NH03. We assess how 
well we have accomplished this goal below. 

We begin by discussing the various bursting regimes and 
the ranges of accretion rates Zacc in which these regimes lie. 
A comparison between the results of ^and Table 3 of NH03 
demonstrates that the two-zone model reproduces all of the 
bursting regimes of NH03 for the relatively high l^^c consid- 
ered in this work, and that the ranges of /acc coincide very 
well with the results of NH03's "canonical" neutron star, for 
which M = 1 AMq, R=10.4 km, and r^oi-e = 10^ K . However, 
the two-zone model cannot reproduce the hydrogen bursts that 
occur for log(Zacc) ^5 ~2.5 because unstable thermonuclear hy- 
drogen burning triggers these burst, but our expression for en 
is temperature-independent. We note that the /acc ranges of 
the various regimes in the two-zone model depend somewhat 
on the values of k and 7He, assumed constant in our model. 

The methods used to determine the equilibrium solution in 
each model differ somewhat, although this difference should 
not affect the eigenvalue. NH03 started with an assumed col- 
umn depth of accreted plasma Eiayei and determined the equi- 
librium solution of the physical quantities p, T, F, X, and Y 
as functions of E. Next, they carried out a linear perturbation 
analysis to determine the stability of the equilibrium solution. 
If 5R(A)facc > 3, where in this case face = S|ayer/S, then the 
configuration was unstable and a type I X-ray burst was as- 
sumed to occur If 5R(A)facc < 3, the configuration was stable, 
and the process was repeated for a larger value of Eiayer- If 
^Wtucc < 3 for all trial values of Siayer, then the system was 
stable, and no type I X-ray bursts occurred. It is important 
to understand that the method of NH03 will find all systems 
that have positive eigenvalues to be unstable because face in- 
creases as Eiayer increases. In the two-zone model, we assume 
complete burning of the accreted matter, and we say that type 
I X-ray bursts occur if 3?(A)facc > and that no bursts occur 
otherwise. If a system is stable to bursts, the accreted mat- 
ter steadily burns to completion, and both methods agree that 
this equilibrium is stable. If a system is unstable to bursts, the 
two-zone model determines that the complete steady-burning 
equilibrium solution is unstable. The method of NH03 may 
find that the identical system is unstable either before all of 
the fuel burns to completion (in a prompt burst, for example) 
or after all the fuel burns to completion (in, say, a delayed 



burst), but it will always determine that the system is indeed 
unstable. Therefore, although the two stability analyses differ 
to some extent, these differences are unimportant with respect 
to the determination of the eigenvalue of a system. 

The eigenvalues derived using the two models agree quite 
well. We have shown already that the real parts of the eigen- 
values agree quite well, for they delineate between the prompt 
and delayed bursting regimes. We now compare the imagi- 
nary parts. We have shown in ^that 3(A) « uju in the de- 
layed mixed burst regime, which means that oscillations oc- 
cur with a period equal to the hydrogen-burning timescale. 
To explicitly show that 3(A) « luu in the global linear sta- 
bility analysis, we consider the calculation presented in §4.2 
of NH03. In that calculation, NH03 consider a system with 
log(/acc) = -0.9, which impUes that E = 1.04 x lO"* gcm"^ s~\ 
since Endd = 8.26 x lO** gcm~^s~' for that system. Figures 
8 and 9 of NH03 imply that E^'' ~ 1-0 x 10** gcm'^, and so 
o^H = 6.5 x lO"'^ s"' from equations i49i and ( 15 0> . This value 
of LUii is very close to the 6.6 x lO"** s"' NH03 quote for the 
imaginary part of the eigenvalue they obtained in their calcu- 
lation. 

Rev nivtsev et al.l J2001I) discovered a class of low- 
frequency oscillations that precede type I X-ray bursts in the 
systems 4U 1608-522, 4U 1636-536, and Aql X-1. These 
oscillations have periods ^120 seconds and have been ob- 
served when the accretion rate Zacc ~ 0.1. While the modes 
resulting from the global stability analysis of NH03 are com- 
plex in this regime, the oscillation periods of these modes are 
roughly a factor of 10 larger than the observed periods. It is 
now clear why NH03 were unsuccessful in explaining these 
observations. We have shown in ^that the oscillation pe- 
riods of overstable modes in this regime are approximately 
equal to the hydrogen-burning timescale tu- From equations 
^ and (EO), 

fH= >770('^V- (54) 

7h(Z^'' + Z„) \Q.7j 

Thus the oscillations preceding a delayed mixed burst have 
too long of a period and cannot be associated with the oscilla- 
tions observed by Revnivtsev et al. (2001) if hydrogen burns 
predominantly via the hot CNO cycle. 

We have shown in §4 that, in the delayed mixed burst 
regime, a considerable amount of stable burning occurs prior 
to a type I X-ray burst. Consequently, a, the ratio of the grav- 
itational energy released via accretion to the nuclear energy 
released during a burst, rises dramatically near the critical 
Zacc above which bursts do not occur. A comparison between 
Figure of this work and Figure 15 of NH03 demonstrates 
that the a-values derived using the two different models agree 
quite well. 

Using an improved version of the model of NH03, 
ICooper et al.l ( 12006.) found that increasing Zq, the CNO abun- 
dance of the accreted plasma, diminishes the range of /acc over 
which delayed mixed bursts occur. In fact, for large enough 
values of Zq, Cooper et al. ( 2006) found that the ranges of /acc 
in which prompt mixed bursts and delayed mixed bursts oc- 
cur are separated by a regime of stable burning. Figure |8j? 
shows the normalized eigenvalues of the two-zone model for 
systems in which Zq > 0.02. The critical /acc above which 
bursts do not occu r decreases with i ncreasing Zq, in agreement 
with Figure 5 of Cooper et al. (2006), but the two-zone model 
does not reproduce the regime of stable burning that separates 
the prompt mixed bursts from the delayed mixed bursts for 
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our choices of k and 7He- However, we are able to repro- 
duce this regime with certain combinations of Zo, k, and 7He- 
Figure Isj? shows the normalized eigenvalues for a model in 
which Zo = 0.132 and 7He = 0.3. A stable burning regime sep- 
arating prompt mixed bursts from the delayed mixed bursts is 
clearly evident for this calculation, showing that the two-zone 
model has the right qualitative behavior. However, our choice 
of Zo for this calculation is more than twice that required in 
the more accurate model of Cooper et al. (2006), showing that 
(not surprisingly) there are quantitative discrepancies. 

5.2. Comparison to Previous Theoretical Work 

For accretion rates /a^c ^0.1, type I X-ray bursts are well- 
understood theoretically, and the results of theoretical mod- 
els are generally in accord with observations. This is not 
true for Zacc ^0.1. Nearly all theoretical models predict that 
accreting neutro n stars should exhibit type I X-ray bursts 
for all /acc < 1 ('Fuiimotoet aL'TosT; "Avasli & Joss" 1981 
Faam 1985; Taam et al. 1996; Bildsten 1998; Fiskeretal. 
12003; Hegeretal. 2005), whereas the model of NH03 and 
the two-zone model presented in this work suggest that bursts 



should cease for 4cc ^, 0.25, 



?reement with observa- 



tions (van Paradiis e Talill979( 1198^ ICornelisse et aljr2003t 
iRemillard et al. 

2006^" Furthermore, our model and that of 
NH03 are the only models able to explain the dramatic in- 
crease in a with increasing /acc observed by v an Paradiis et al.l 
il988h . In this section, we compare and contrast the two-zone 
model presented in this work with previous one-zone burst 
models. Additionally, we discuss the results of this work in 
relation to those of detailed time-dependent multi-zone mod- 
els that utilize sophisticated nuclear reaction networks. 

We begin by discussing the stabihty criterion of one-zone 
models (e.g., Fuiimoto et al. 1981; 'Pacz vnski' 1 98 3 a^, 'Bildsten] 
[1998; Cummins & Bildsten 2000; Heger et al. 2005). These 
models presume that helium burning (or hydrogen burning for 
very low Zacc) initiates a thermonuclear instability and pro- 
duces a type I X-ray burst. As mentioned in ^ the hydro- 
gen burning rate is temperature-independent for T > 10^ K 
and so hydrogen burning cannot be thermally unstable for 
hcc ^ 0.01. Since all of the one-zone models give very similar 
results, we focus on the model of Heaer et al. (2005), which is 
an extension of the model of Paczvhski ( 1983a|), for a compar- 
ison. Written in our notation, the differential equations (1-2) 
of Heger et al. ( 2005) that govern their model are 



dt 

dt 



eHe(7He) 



Co 



eHe(7He) 



'He 



(55) 



(56) 



where = acT^^/liKYiii^. They presume that Jf^ > S^'j. and 
that no stable helium burning occ urs pri or to the burst. The 
equilibrium solution to equations j55l56> is 



eHe(7He)SHe - ^He - ^^^HeJ 



or equivalently. 



eHe(7He) - Ecool, 



(57) 
(58) 

where EcooI is given by equation (15 1> . Using the thin-shell 
thermal instability criterion of equation i52i . one finds that 
the transition from stability to instability occurs when ly = 
91neHe/91nr w 4, since 91necooi/91nr = 4. One finds nearly 
the same condition by performing a standard linear stability 
analysis on equations (I55> and (I56> . Thus, from equation 



il5\ . the equilibria are stable to bursts when T^'^ 6 x 10** 
K. However, Figures|2land|3illustrate that many of the equi- 
libria in our two-zone model are stable to thermal perturba- 
tions for temperatures much less than 6 x 10^ K. To under- 
stand this, note that equation ( I58> implies that helium burn- 
ing is the only source of nuclear energy generation. In our 
two-zone model, both hydrogen burning and helium burning 
are sources of nuclear energy generation, and so qualitatively 
ficooi ^ ^HeiT) + euiZ). Hydrogen burning is temperature- 
independent, so equation i52\ still gives the right criterion for 
a thin-shell thermal instability. However, because Ccooi > fHe, 
the critical value of i/ for instability is larger than 4. Roughly, 
on expects something like i^cnt ^ 4(ecooi/eHe)- Several authors 
(e.g., Fuiimoto et al. 1981; Cummins & Bildsten 2000) have 
included the effect of hydrogen burning on the stability cri- 
terion in their one-zone models. These models assume that 
no helium burning takes place prior to ignition, so Z = Zo 
throughout the accreted layer. Consequently, they find that the 
effect of hydrogen burning on the stability of the system is mi- 
nor In contrast, stable helium burning prior to ignition is sub- 
stantial for /acc ^0.1 in both our two-zone model and that of 
NH03, and so the effect of hydrogen burning on the stability 
of the system is considerable. Thus, hydrogen burning via the 
temperature-independent hot CNO cycle, augmented by the 
extra CNO produced from stable helium burning, helps sta- 
bilize nuclear burning on accreting neutron stars even at tem- 
peratures 3 X 10** K, which is much less than the 6 x 10** 

K needed to stabilize pure helium burning. 

In the limit Xo and Zo 0, equations ( I31I35> of our 
two-zone model reduce to the following set of two coupled 
differential equations: 



dY. 



He 



dt 



eHe(7He) 



dt 



4Co 



eHe(7He)" 



He- 



''He 



'He 



dH 



He 



4I]He dt 



(59) 



(60) 



where again Fue = acT^^/3KY.ue from equation and the 
equilibrium is given by equation i51\ . Equations ( I59l l and 
i60\ are effectively the sa me as eq uations ( I55t and ( I56t . With 
the system of equations ( I59I60> . we are able to reproduce 
the results of the o ne-zone models of lPaczvri ski ( 1983a) and 
iHeger et'al] J2005ft . As expected, the critical temperature for 
stabihty is ~ 6 x 10^ K in this case. 

It is perhaps not surprising that one-zone models that fo- 
cus only on helium burning are not able to reproduce all of 
the bursting regimes found using global linear stability anal- 
yses or observed in nature. As we have shown in this pa- 
per, one needs at least two zones, one in which only he- 
lium burns and another in which both hydrogen and he- 
lium burn together, to incorporate all of the key physics in- 
cluded in the model of NH03. One-zone models simply 
do not include enough physics. More troubling is the fact 
that the results of d etailed multi-zone calculations of type 
I X-ray bursts (e.g.. 'Avasli & .Toss"1982t iTaam et alJ 11995 
Fiskeretal. 2003; Heseretal. 2005) do not reproduce the 
regime of delayed mixed bursts and are also inconsistent with 
observations for Zacc ?i 0.1. Since we do not have access to 
such models, we are unable to conduct a direct comparison 
between our model and these numerical calculations. How- 
ever, we attempt to better unde rstand this discrepancy by con- 
sidering the work of Fisker et al. (2006). 

Fisker et al. (2006) study the effects of the relatively un- 
constrained '^0(a,7)''Ne reaction rate on type I X-ray 
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bursts from neutron stars accreting at 4(.c ~ 0.1. This re- 
action is one of the two hot CNO cy cle breakout reactions 
JWagoneJ fT969l IWallace & Wooslevlf!9 81: Langan keetall 
ll986H Wiescher et al.lll999 ^.'Schatz et aDllQQQ ). For the low- 
est of the three trial reaction rates they consider, they find 
that, after an initial type 1 X-ray burst, hydrogen and helium 
burn stably via the hot CNO cycle and triple-a reactions, re- 
spectively, and the nuclear burning generates a slowly oscil- 
lating luminosity with a period approximately equal to the 
hydrogen-burning timescale. The nuclear burning behavior 
in this calculation is nearly identical to that of the nuclear 
burning preceding a delayed mixed burst found in ^ This 
similarity is perhaps to be expected, for the limited reaction 
networks employed in the model of NH03 and the present 
model omit hot CNO cycle breakout reactions en tirely. For 
the two higher trial reaction rates investigated by Fisk er et al] 
( 12006) . the breakout sequence '^0(a,7)'^Ne(/?,7)^"Na dimin- 
ishes the CNO abundance and thus reduces the rate of hydro- 
gen burning. Hence, when helium ignition commences, the 
helium energy generation rate dominates the total nuclear en- 
ergy generation rate, and a prompt type 1 X-ray burst occurs 
via a thin-shell thermal instability. We tentatively suggest that 
the discrepancies between the results of the global linear sta- 
bility analysis of NH03 and the results of multi-zone calcula- 
tions arise from differences in the treatment of hot CNO cy- 
cle breakout reactions. In the NH03 model and our two-zone 
model, hydrogen burning via the temperature-independent hot 
CNO cycle helps stabilize nuclear burning on accreting neu- 
tron stars. Breakout reactions would reduce the degree to 
which hydrogen burns via the hot CNO cycle and thereby in- 
crease the temperature-sensitivity of the total effective nuclear 
energy generation rate. Thus, breakout reactions would pre- 
sumably extend the range of accretion rates in which nuclear 
burning is thermally unstable, but we do not include these re- 
actions. That observations agree much better with the results 
of NH03 and our two-zone model may imply that the true 
cross sections of reactions such as the hot CNO cycle breakout 
reactions '^0(Q;,7)'^Ne and '^Ne(a,/?)^'Na are much smaller 
than the cross sections employed in the reaction networks of 
multi-zone models, i.e., the true '^0(Q;,7)'''Ne rea ction rate is 
closer to the lowest reaction rate considered by Fisk er et al] 
(I2OO6). By omitting hot CNO cycle breakout reactions alto- 
gether, the model of NH03 is perhaps a better representation 
of the nuclear burning that precedes type 1 X-ray bursts than 



time-dependent multi-zone models as presently implemented. 

The model we present in ^is valid only when Sh < Shc, 
which restricts the range of accretion rates we can study to 
Zacc ^0.3. We have constructed a similar model to deter- 
mine the stability of nuclear burning on accreting neutron 
stars for which Eh > Shb, allowing us to study the stabil- 
ity of nuclear burning at higher l^^c- Although we find that 
all equilibria for Zacc ^ 0.3 are stable, which is consistent with 
the results of NH03, we cannot state with confidence that ei- 
ther the two-zone model or the model of NH03 is an accu- 
rate representation of nuclear burning on accreting neutron 
stars for such high accretion rates. For Zacc 1, nuclear re- 
actions other than the hot CNO cycle and triple-a reaction 
almost certainly play a significant role in the nuclear burning 
prior to a type 1 X-ray burst. It is possible that there exists an- 
other unstable burning regime near l^^c ^ 1 , and the range of 
Zacc in which this regime might exist would probably depend 
upon the /acc at which the effective hydrogen-burning rate 
becomes predominantly temperature-dependent rather than 
composition-dependent. While most low-mass X-ray bina- 
ries with Zacc ^ 0.3 do not have type 1 X-ray bursts, GX 17+2 
and Cyg X-2 are notable exceptions cKahn & Grindlav|l[1984t 
Tawara et a l]|1984t lSztaino et alJI 198 61: iKuulkers et alJI 19951 
1997; WiinandsetalJ 119971 iSmald 119981 IKuulkers et alJ 
2002), although these sources may exhibit type I X-ray bursts 
for other reasons. For example, bursts could occ ur at l„c.c ^ 1 if 
the ac creted plasma is hydrogen-deficient (e.g. JCooper et alJ 

6. CONCLUSIONS 

We have constructed a simple two-zone model of type 1 X- 
ray bursts on accreting neutron stars. This model reproduces 
the delayed mixed burst regime of NH03 as well as the helium 
and prompt mixed burst regimes of previous studies (Fuji- 
moto et al. 1981; Fushiki & Lamb 1987; Cumming & Bildsten 
2000; NH03) , and it agrees well with observations of type I 
X-ray bursts (van Paradiis et al. 1979'. 1988; ' Cornelisse et alJ 
^03; Remillard et al. 2006). More importantly, the model il- 
lustrates the physics of the onset of instability as a function 
of the local accretion rate E, and it facilitates comparisons 
between global linear stability analyses (which are more ac- 
curate but difficult to understand physically) and other burst 
models. 

A pure, rapidly growing thermal instability in the heUum- 
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burning zone triggers bursts at relatively low S. As S in- 
creases above O.lSEdd, the trigger mechanism evolves from 
that of a thermal instability to that of a slowly growing over- 
stability involving all parameters, particularly the CNO mass 
fraction Zh of the hydrogen-burning zone. The competition 
between nuclear heating via the /3-limited CNO cycle as well 
as the triple-a process and radiative cooling via outward dif- 
fusion of photons coupled with radiation from the stellar sur- 
face drives oscillations with a period approximately equal to 
the hydrogen-burning timescale /E. If these oscillations 
grow in time, the temperature Tm at the base of the helium 
layer will rise, eventually triggering a thin-shell thermal insta- 
bility and hence a delayed mixed burst. For S/SEdd ^ 0.25 
radiative cooling from the stellar surface dampens the oversta- 
bility, and no bursts occur. We consider our two-zone model 
to be "minimal," by which we mean that no model consisting 
of only a proper subset of the five time-dependent variables 
Eh, Zh, 7h, Ehb, and Tne is able to reproduce the phenomenon 
of delayed mixed bursts. 

Nearly all other theoretical models predict that bursts 
should occur for all E/EEdd 1^ in disagreement with the 
results of both NH03 and the two-zone model, as well as with 
observations. We suggest that this discrepancy arises from 
the assumed strength of the hot CNO cycle breakout reac- 
tion '^0(a,7 )''^Ne (Wagoner 1969; Wallace & Woosley 1981; 
ISchatz et al.l [T999t iFisker et al.i ,2006) in time-dependent 
multi-zone burst models. That observations agree much bet- 
ter with the results of both NH03 and this work may im- 



ply that the true '^0(a,7)'''Ne cross section is much smaller 
than the cross sections employed in the reaction networks of 
these models. Further calculations such as those presented by 
Fisker et al. ( 2006) in which the '^0(a,7)'''Ne reaction rate in 
the networks of time-dependent multi-zone models is varied 
should be performed. 

We have considered only two forms of nuclear burning in 
our model: hydrogen burning via the hot CNO cycle and 
helium burning via the triple-a process. While this sim- 
plification is almost surely reasonable for studying the on- 
set of helium bursts, certain additional reactions may be im- 
portant in the delayed mixed burst regime or at yet higher 
E. In particular, hot CNO cycle breakout reactions such as 
'^0(a,7)'''Ne could significantly affect the CNO metallicity 
Zh of the hydrogen burning zone and consequently have some 
effect on delayed mixed bursts. Additional hydrogen burn- 
ing processes that could circumvent the hot CNO cycle such 
as the i50(p,7)(/3+zy)i''0(p,7)i^F(p,7)i'*Ne(/3+i/)i'^F(p,a)i50 
reaction sequence (Stefan etal. 2006) could affect delayed 
mixed bursts as well. These issues are worthy of further in- 
vestigation. 
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